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ABSTRACT 



A coil antenna consisting of approximately 6000 turns 
of copper wire was utilized to measure the horizontal com- 
ponent of fluctuations of the earth's magnetic field on the 
floor of Monterey Bay in water depth of approximately 50 
meters. The results indicate that the power spectral den- 
sity of the fluctuations varies from lOnT^/Hz at 0.04 Hz to 
10“^nT^/Hz at 25 Hz, a monotonic decrease of about 6 dB/ 
octave, except in the 8-20 Hz region where the Schumann 
resonances occur. While the sensitivity of the equipment 
was insufficient to measure the vertical component of the 
fluctuation we can put an upper limit of 10“3nT2/Hz at 1 Hz 
and 10“^(nT)2/Hz at 10 Hz on the magnitude of this component. 

The same sensor was also used to measure various com- 
ponents of the field fluctuations at a remote land' site 
(Chew's Ridge). In the frequency range observed the gen- 
eral shape of the spectra was similar to those obtained at 
sea. However, a strong azimuthal variation at certain 
frequencies was noted in the land data. The possibility 
that these directional signals are of man made origin 
cannot be excluded at this time. 
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I. INTRODUCTION 



This thesis research is a part of an ongoing effort i at 
the Naval Postgraduate School, to obtain improved long term 
data and interpretations of the electromagnetic noise on 
the ocean floor. The overall project objectives are to 
study and interpret signals covering 5 decades of frequency 
from .001 Hz to 100 Hz through the combined use of total 
field magnetometers and ULF/ELF receivers. The overall pro- 
ject emphasizes the importance of obtaining measurements of 
geomagnetic noise fluctuations on the sea floor and on land 
over a period of several years and at various locations. 

The particular objectives of this thesis research are 
1) to collect local magnetic field fluctuation data at a 
relatively quiet land site near Monterey, California for 
comparison with data obtained from the floor of Monterey 
Bay and 2) to test a coil type antenna specifically de- 
signed for deep sea use, as well as the required electronics 
and data recording equipment. 
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II. BACKGROUND 



A. SOURCES OF GEOMAGNETIC FLUCTUATIONS 

Geomagnetic fluctuations result from sources internal 
and external to the earth's surface. Internal sources, 
which account for approximately one-third of the variations 
in earth's magnetic field, are referred to as secular vari- 
ations, having periods of thousands of years. Geomagnetic 
storms, diurnal variations, and micropulsations are due to 
external sources and have periods on the order of days and 
less. 

Probably the best known sources of geomagnetic fluctua- 
tions are those generated by the interaction of charged 
particles with the magnetic fields of the magnetosphere. 
Figure 1 indicates an approximate power spectrum of vari- 
ations measured on the earth's surface having periods of a 
day or shorter. Although this phenomenon is not fully 
understood, much effort has been expended in studying the 
sources and predicting the occurrences of these variations. 
The high frequency fluctuations, referred to as micropul- 
sations, fall into two general classes: Pi or irregular 
pulsations and Pc or continuous pulsations and are be- 
lieved to result from wave- particle interaction in the 
magnetosphere (Jacobs 1970). 

Of particular interest are long and short term predic- 
tions of the Pc 1 geomagnetic pulsations (.2-5 Hz). A 
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POWER DEl'JSITY ( nT ) 




Frequency (Hz) 

Figure i 

Power Spectrum of Geomognetic Disturbances Observed 
on the Surface of the Earth. 

[ Clad is, Davidson, and Newkirk 1971 ] 
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number of methods are available for the forecasting of long- 
term (Fraser-Smith , 1972 and Matyeyeva et al . , 1972) and 
short-term (Fraser-Smith, 1980) geomagnetic pulsation 
occurances. However, the accuracy of these methods must 
eventually be determined by comparison with the observed 
data for both research and operational reasons. 

Thus a better understanding of the time dependent 
component of the geomagnetic field is of considerable 
interest both in geophysical and naval applications. Pre- 
diction of geomagnetic fluctuations is especially valuable 
in Navy operational and projected systems such as extreme- 
ly low frequency (ELF) communications, submarine detection 
(Magnetic Anomaly Detection), and weapon deployment (mines, 
particle beams, and torpedoes). 

In addition to the information about magnetospheric 
processes which can be deduced from a study of time varia- 
tions of the geomagnetic field, these fluctuationa also 
serve as an important tool in the study of other geophysical 
quantities. Thus both conventional and naval magnetotelluric 
techniques make use of the variable geomagnetic field in 
the measurement of ocean floor conductivities. Similar 
measurements are employed in the study of earthquake faults, 
the location of aquifers and geophysical prospecting in 
general . 
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B. ELECTROMAGNETIC PROPAGATION IN LAYERED CONDUCTING MEDIA 
To describe electromagnetic propagation in layered con- 
ducting media, a three layer model was selected (see figure 
2) . Note that plane wave propagation is assumed. This mod- 
el was used to predict the expected geomagnetic variations 
in the media and at their interfaces . Before proceeding 
with an explanation of the model and its use, the follow- 
ing definitions are presented: 

X = Wavelength, meter 

V = Phase velocity, meter/second 
P 

r = a+j3. Complex propagation constant, meter 
CO = Angular frequency, radian/second 
U = Permeability, henry /meter 

e = Dielectric constant, farad/meter, e=£ £- where £«= 

_ 9 r u u 

10 /36it farad/meter 

£ = Relative dielectric constant, for air £ =1 and for 
r r 

sea water £ =81 
r 

a = Conductivity, siemens/meter, for sea water <^2=4 

siemens/meter and for sea floor siemens/meter 

E = Electric field strength, volt/meter 

B = Magnetic flux density, tesla 

6 = Skin depth, 1/a, meter 

f = Frequency, hertz 

n = Intrinsic impedance, ohm 

The media are assumed to be non-magnetic with V* 2 =^ 12 = 

_7 

li2=ViQ=4TT (10 ) henry /meter. Sea water and the ocean floor 
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are assumed to be good conductors in the frequency of in- 
terest (.04 to 25 Hz), i.e., a/coe>>l (conduction current 
dominates) . On the other hand, air is assumed to be a good 
insulator, i.e., a/coe<<l (displacement current dominates). 
For these conditions the following equations for the propa- 
gation constants and intrinsic impedances along with their 
numerical values are valid: 

Region 1: (air) 

r^=jo) (y^e^)^ I r [ =1. 9 (10~^) f meter ^ 





1=377 ohm 




Region 2: (sea 


. water) 






, h 1 




3v 


-1 


r2=(ja)yQcr2) | 


121=5.6(10 


J) f 2 


meter 


Tl2= 1 


ri2l=1.4(10” 




ohm 


Region 3; (sea 


, floor) 






% . 


1 “ 


4 h 


-1 


r3=(3wyQa3) | 


131 = 8 . 9(10 


)f ^ 


meter 


n3=(jwyoA3)^ i 


7131 = 8 . 9 ( 10 " 


^)f^ 


ohm 



For simplicity, only the theory for a normally inci- 
dent, plane polarized wave is considered. The modification 
for the general case of oblique incidence is easily made. 

It can be shown that the direction of propagation in the 
sea is vertically downward, regardless of the angle of 
incidence . 

Use of Maxwell's equations and the above assumptions 
results in the following wave equations for the media of 
interest; 
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Figure 2. THREE-LAYER MODEL 



Region 1: (air) 

E = E.e“^^^+E 
XI 1 r 



B 



'yi= {Uq/ti^) 

Region 2: (sea water) 

-E 



JlZ- 



E = E.e‘^2^-- 



X2 



Region 3: (sea floor) 



E = E e~^®^ 

X3 t 

®y,= <“o/^3>V 



-Taz 



The boundary conditions require that the electric and 
magnetic fields be continuous at each • interface , i.e., at 
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the air-sea water and sea water-sea floor boundaries. 
These conditions are: 



z= 0 : E.+E = E^+E 

1 r + - 






-Fad, 



z=d; E.e " 2 ^+E e^ 2 ^ = E^e 



-Fad 






-Fad 



“t 

2d_^:, ^Fadj ^ 

We have here five unknowns with only four relations 
between them. To obtain a meaningful result from the above 
equations, the various E components are expressed as 
ratios, normalized to E^. The form of the ratios are: 

E^ = (E^E^), E^ = (E_^/E^), ^ = (E_/E^) , and E^ = 
(E^/E.) 

Algebraic calculations result in the following expressions: 

2n 2 (n 2+n 3 ) 

E = - 2 Fad 

— (• (ns+na) (ni+n2) + (n2-ni) (n3-n2)e 

2 n 2 (n 3 -n 2) 



^ (na+ha) (h i+ri 2 ) (h 2“n I ) (n3~'H2) 

Et = E^e-<''^-''=>^E.e"'2+r,)a 

E = E^+E_-l 

Using the above equations for region 2 (sea water) , 
an expression for the induced-to-incident magnetic field 



on the sea floor is obtained. 

B (z) = E /na = E. (^e"^^^-E e^^^) since E . = (n i/Uo ) B . 

(z) = (n i/yo )BjE^e"^ 2 Z_g ^FjZj 

ya 1 — ^ — 

(z) = B /B. = (ni/yo){E^e”^ 2 ^-E_e^"^} 

— y 2 y 2 1 I- 

Now assuming that ri 2 and n 3 are much less than ri i , we define 
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q = Tia/na = { 02 / 03 }^. After some algebraic manipulations, 
we obtain: 

1 s (d-z) , . -T 2 (d-z) 

B_. (z) = 2{i2±l)e llaille j 



-Yz 



/ _Li\ f2d f ~F2d 

(q+l)e -(q-l)e 



For z=d and substituting the numerical values of r\ 2 i Hs/ 
F 2 into the above equation at a specified frequency, one 
can obtain a plot of the induced-to-incident magnetic 
field versus depth. See Figure 3- Clearly, if one is to 
detect magnetic field variations in sea water at any 
appreciable depth, then the variations must be of low 
frequency. 
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RATIO OF INDUCED TO INCIDENT FIELD AMPLITUDE 




Figure 3. 
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C. GEOMAGNETIC FORECASTS 



The Space Environment Services Center located in 
Boulder, Colorado publishes a Preliminary Report and Fore- 
cast of Solar Geophysical Data, including geomagnetic 
activity. A more detailed synopsis of geophysical data is 
contained in the environmental data services monthly publi- 
cation "Solar-Geophysical Data". Of importance to the re- 
search effort here at Naval Postgraduate School are the 
reports of the A and K-indices for the mid-latitude region 
containing Monterey. The major monitoring station for the 
mid-latitude areas in the U.S. is located in Fredericks- 
burg, Virginia. 

A and K-indices reported by Fredericksburg are of in- 
terest to NPGS, since this station lies approximately on 
the same latitude as Monterey. The A-indices indicate a 
24-hr daily variation in the earth's magnetic field while 
the K-indices indicate a 3-hr period of variation within 
the 2^-hr A-index. The K-indices vary in range from 0 
(very quiet) to 9 (extremely disturbed) with an A-index of 
30 or greater indicating geomagnetic storm conditions. 

The above information was utilized to correlate geo- 
magnetic predictions and forecasts with geomagnetic data 
obtained in Monterey Bay. This correlation concept is util- 
ized below in Section V, experimental results. (See Appen- 
dix C for a further description of geomagnetic forecasts). 
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D. REVIEW OF EARLIER WORK 



A fair amount of research has been and is being done in 
the detection and prediction of low frequency variations 
(less than 100 Hz) of the earth's magnetic field, both on 
land and in the sea. Earlier experiments by Horton and 
Hoffman (1962), Santirocco and Parker (1963)> and Davidson 
(1964) revealed that the power spectral density of back- 
ground geomagnetic activity in the frequency range 10“^- 1 
Hz decreased monotonically with a slope of approximately -6 
dS/octave. The spectrum then levels off and is dominated 
by the Schumann resonances in the frequency range l-40 Hz 
(Schumann and Konig, 1954) <. Measurements indicate a gen- 
eral increase in background activity between 20 and 70 Hz 
as reported by Campbell (1966), with a continued sharp 
decrease in slope above 80 Hz. 

Spectrums obtained by Larsen and Egeland (1968) indi- 
cate that the general decline observed at frequencies below 
5 Hz continues through 70 Hz and probably to higher fre- 
quencies (Maxwell and Stone, 1963; Westerlund, 1970) with 
the Schumann resonance activity being superimposed on the 
declining background as a form of fine structure. Hence, 
there is a remarkably steady decline of the background fluc- 
truations at a rate of approximately 6 dB/octave covering 
almost 6 decades of frequency (10“^ - 10^ Hz) (Fraser-Smith 
and Buxton, 1975) » 
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Recent measurements conducted by Fraser-Smith and Bux- 
ton in the .1 - Hz frequency range confirmed earlier 
predictions and detections. They discovered that the gen- 
eral decline in background activity continued for frequen- 
cies up to 5 Hz , but for frequencies above 5 Hz the decline 
appeared to be suppressed by the occurence of the Schumann 
resonances. 

More recent work was accomplished in the low frequency 
ranges at La Mesa in Monterey. Barry (1978) and Clayton 
(1979) investigated geomagnetic fluctuations in the fre- 
quency ranges .1-10 Hz and .4-40 Hz, respectively. Their 
results are in general agreement with earlier measurements. 

Chaffee (1979) measured geomagnetic variations of the 
earth's magnetic field on the floor of Monterey Bay in 100- 
200 meters of sea water and covering the frequency range 
,01-3 Hz. Chaffee observed the power spectrum of typical 
geomagnetic fluctuations to vary from 1 nT^/Hz at ,01 Hz to 
approximately 10“5 nT^/Hz, Most low frequency studies made 
in sea water have been performed to detect geomagnetic 
fluctuations produced by surface waves. Since little data 
is available on actual geomagnetic noise at low frequen- 
cies, the project at Naval Postgraduate School is a con- 
tinuing effort to create a reliable data base for future 
reference and correlation. 
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Ill . EXPERIMENTAL EQUIPMENT AND TESTS 



A. EQUIPMENT CONFIGURATION 

Most of the equipment utilized for this project was 
either designed and/or manufactured at the Naval Postgra- 
duate School, while test equipment used was mostly of 
commercial origin. Much effort was expended in testing and 
ensuring that the experimental equipment operated as con- 
cieved and predicted. It was important for all equipment 
to be dc coupled to avoid loss of low frequency signals. A 
discussion of the equipment and its configuration follows: 

1. Data Acquisition Equipment 
SENSOR 

The general hookup of the data acquisition 
equipment is depicted in Figure 4. The sensor is a self- 
supporting, continuously wound, non-center-tapped coil an- 
tenna manufactured from approximately 6000 turns of 18 gauge 
copper magnet wire. The sensitivity of the coil is estab- 
lished in Section III. B. Its weight is approximately 100 
pounds. Its overall dimensions (.15ni tiy .37m) were deter- 
mined by the largest glass sphere commercially manufactured 
by Benthos, Incorporated. Refer to Figure 5* The coil re- 
sistance is 120 ohms and its self-inductance is approxi- 
mately 9*31 henries. 
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SENSOR SPHERE INSTRUMENTATION SPHERE 



(LOVJER) (UPPER) 

1 I 1 




I I 



Figure 4 Data Taking Equipment 







Figure 5 



Sensor Dimensions 



PREAMPLIFIER 



Various differential preamplifiers were tested 
and used including the PAR-113 amplifier manufactured by 
Princeton Applied Research Corporation. This preamplifier 
was selected because of its higher signal to noise ratio, 
low power requirements r and compactness. The preamplifier 
itself is packaged as an integral part of the sensor, i.e., 
it is housed in the glass sphere with the sensor. Attempts 
to separate the preamplifier from the sensor resulted in 
additional noise caused by cable motion. 

The front end of the preamplifier utilized a low- 
noise, "state of the art" differential device (LH OO 38 C) 
manufactured by National Semiconductor Corporation. It has 
a single-ended output which is connected to additional 
stages of amplification. The final stage of the preampli- 
fier contains an active filter which is set to a cutoff 
frequency of 20 Hz. The overall preamplifier gain, in- 
cluding filters, is approximately 66 dB. A graph of the 
preamplifier gain characteristics is included as Figure 14. 

A schematic diagram of the preamplifier circuit is contained 
in Appendix A. 

VOLTAGE CONTROLLED OSCILLATOR 

The schematic diagram of the voltage controlled 
oscillator (VCO) is included in Appendix A. The function 
of the VCO is to convert low frequency analog voltages re- 
ceived from the sensor through the preamplifier into a 
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varying frequency, centered at I 5 OO Hz. The output level of 
the VCO is approximately 0.2 volt (zero to peak) square wave 
of varying frequency which is recorded on an analog cassette 
tape recorder. The conversion factor for the VCO is 1.25 
Hz/mV. ■ 

REFERENCE OSCILLATOR 

A 2 KHz reference oscillator signal is recorded 
simultaneously with the VCO signal on a separate channel. 

The schematic diagram of the reference oscillator is also 
included in Appendix A. The primary purpose of the refer- 
ence oscillator is to reduce tape recorder noise (wow and 
flutter) during playback. 

TAPE RECORDER 

The tape recorder employed for this project is a 
portable stereo cassette type (Marantz Superscope, Model 
CD- 330 ) » powered by a 6 Vdc battery. The Superscope was 
selected to meet the size constraints imposed by the glass 
spheres. Additional characteristics are: 

(1) Tape Speed: l-?/8 in/s 

( 2 ) Frequency response: 40 Hz to 12 KHz 

( 3 ) Signal -to -noise ratio; 50 dB 

(4) Wow and flutter: .12^ 

POWER SUPPLIES 

Gould Gelyte lead batteries are utilized to meet 
dc power requirements. They are sealed, rechargeable, and 
virtually maintainance free batteries, containing no free 
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electrolyte. Various sizes are used to meet the power 
specifications of all data taking and analysis equipment, 
except test equipment . 

2 . Data Analysis Equipment 

Refer to Figure 6 for the general hookup of the data 
analysis equipment. The tape recorder and power supplies 
(discussed above) were employed in data analysis as well as 
in data taking. A description of the additional analysis 
equipment follows: 

MIXER 

Appendix A contains the schematic for the mixer 
circuit. The mixer is powered by t 12 Vdc and accepts the 
recorded sensor data and reference oscillator signals at its 




Figure 6 DATA ANALYSIS EQUIPMENT 

inputs from the tape recorder. The difference frequency of 
the mixer is the input to a frequency-to -voltage converter. 
The mixer circuit also contains a low pass filter with a 
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cutoff frequency at approximately 1200 Hz. The purpose of 
the low pass filter is to elimate unwanted harmonics of the 
difference frequency. 

FREQUENCY -TO -VOLTAGE CONVERTER (ANADEX PI-375) 

The purpose of the frequency-to-voltage converter 
(FVC) is to convert the difference frequency from the mixer 
into an analog voltage. The varying voltage is an exact 
representation of the geomagnetic fluctuations detected hy 
the sensor. The FVC is powered by 2k Vdc and has the 
following additional characteristics: 

(a) Input frequency range: 0-2000 Hz 

(b) Output dc voltage: O-IOV full scale 

(c) Conversion factor: 4.91 mV/Hz 

The output of the converter is analyzed with a spectrum 
analyzer. 

SPECTRUM ANALYZER (S/A) 

Numerous spectrum analyzers were used for data 
analysis: The Schlumberger Model 1510-03» The Ubiquitos 

Model UA-500, The Mini-Ubiquitos Model 440, and the Hewlett 
Packard Model 3582A. Each of these analyzers has inherent 
advantages and disadvantages revealed during their operation. 

The Schlumberger S/A does not provide for dc 
coupling and hence was not adequate for low frequency appli- 
cations. The Ubiquitos has dc coupling and limited low 
frequency capabilities. Also an external display had to be 
provided to obtain a power spectral density read out which 
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would be highly inaccurate. The Hewlett Packard' Analyzer 
was capable of dc coupling and provided for a digital dis- 
play of the power spectrum. However, it does not have a 
logarithmic scale for frequency. 

The most adequate S/A utilized was the Mini- 
Ubiquitos Model 440. It provides for dc coupling, has a 
digital display of the power spectrum, and has a logarithmic 
as well as a linear frequency scale. 

Since the discriminator output is a positive dc 
voltage, it was necessary to provide for dc offset in order 
to increase the input sensitivity of the spectrum analyzer. 
The Tertronix Model AM ^02 differential amplifier was 
utilized with its gain set to XI. 

3 . Underwater Equipment 

Underwater equipment included the following (see 
Figure 7 ) '• 

(a) Benthos vacusealed glass spheres v;ith 
associated hard hats, penetrators, and sealant strip. 

(b) Cable and connectors 

(c) Anchor and nylon line 

The particular spheres employed during this project were a 
low expansion borosilicate type, having a 0.404 meter in- 
side diameter. The actual spheres used in the tests were 
two hemispherical shells specially matched to guard against 
leakage. They are pressure tested by the manufacturer, to 
1000 atmospheres, which is equivalent to a sea water depth 
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SCHEMATIC OF EXPERIMENTAL LAYOUT 



FOR AT SEA MEASUREMENTS 




FIGURE 7 
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of 10,000 meters. Before the spheres are placed in the 
hardhats, which provide protection against breakage, the 
hemispherical seam is sealed with sealant strip and taped. 

Penetrators were installed in the wall of the 
Benthos Glass Housings to provide electrical connections to 
enclosed equipment. A pressure type penetrator is used to 
activate the tape recorder when the instrumentation sphere 
is at a depth of approximately 10 meters. 

Signals from the sensor to the instrumentation 
sphere are propagated through a twisted-pair, shielded, 
rubber coated, 25 meters long cable with Brantner type sea 
connectors. The connectors are made from a special molded 
neoprene material; they have brass, hard gold plated con- 
tacts and are pressure rated to greater than 1000 atmos- 
pheres when mated. 

An anchor supports the sensor sphere and provides 
added negative buoyancy to prevent the system from floating. 
It was designed and constructed from PVC piping, filled 
with sand, and neoprene-coated lead bricks were strapped to 
the piping for additional weight. Extreme care was taken 
to ensure that any exposed metallic surface which might 
induce unwanted currents into the sea water are eliminated. 
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B. SENSOR SENSITIVITY 

1. Theoretical Sensor Sensitivity 

The use of wire coils to measure magnetic field 
changes is not new. Faraday's Law of Induction states that 
the electromotive force (emf) induced in a single-turn 
circuit is equal to minus the rate of change of flux: 

emf ="|^ (3-1) 

The flux is that which passes through a surface 
bounded by the conductor and is defined as: 

$ = /B-da (3-2) 

If we consider N turns of wire connected in series, 
then the total induced emf is the sum of the emf induced in 
each separate turn, or for a loop: 

emf = -N II (3-3) 

In this work it shall be assumed that magnetic 
fields vary sinusoidally with time and may be written as: 

B = B e"^“^ (3-4) 

o 

Where B is a magnetic vector field and only a 
o 

function of spatial coordinates and not of time. By substi- 
tuting equation (3-2) into (3~3) and with equation (3-^) we 
have: 

emf = -NB jwA (3-5) 

o 



33 



where A is the area enclosed by the coil. The magnitude of 
this emf is; 



emf = NB (2Trf)A (3-6) 

o 

If the -output of the sensor is connected to a high input 
impedance amplifier, then negligible current can flow 
through the sensor which might otherwise produce a back emf 
and reduce the magnitude of the field to be measured. The 
sensor's open circuit voltage or sensitivity can be calcu- 
lated from the parameters given in Section III. A. 1 for a 
1 nT field at a frequency of 1 Hz; emf = 3 ‘ill? V'V, thus 20 
log (emf) = -110ol4 dB. (re: 1 Volt) 

The area used in calculating the emf for a single 
turn of wire is easier to visualize than the area of a 
finite width coil. It was decided that an average area 
obtained from the area within the inner radius and the area 
within the outer radius would be used. This method appears 
to be more accurate than using an average radius to compute 
the area. Figure 8 depicts the expected sensitivity of 
the sensor. 
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2. Experimental Determination Of Sensor Sensitivity 
Two methods were utilized to determine the sen- 
sor's actual sensitivity. The first method consisted of 
producing a known field at a specified distance from an 84 
turn copper coil and measuring the electromotive force or 
voltage signal produced by the sensor. The second method 
was similar to the first, but a .61 meter radius Helmholtz 
coil was used to establish a more uniform magnetic field in 
the central region. 

METHOD I. 

The magnetic field produced by a circular ring of N 
turns at a point Z along the symmetry axis is given by: 

(z" + rM 



Equation (3“7) is obtained from the Biot-Savart Law using 
the units of I (ampere), r (meters), Z (meters), and B 
(nanotesla) as depicted in Figure 9» Test coil specifica- 
tions were as follows: 
r = .1577 m 
N = 84 turns 
r inner = 15 cm 
r outer = 16. 5 cm 

where r was computed from the average area of the coil. 



Rearranging equation (3~7) we have; 

B. 



(z^ + r") 
Irms = 



2 \ ^ / 2 



(3-8) 



1.68 (10‘*)r2TT 
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Specifying the desired field strength and distance deter- 
mines the current through the coil at different frequencies. 
To produce the required current a llOfi resistor was placed 
in series with the coil of depicted in Figure 10, and 
the signal voltage was set according to: 




I(R + R ) 
1 2 



(3-9) 



The current was determined to he within 1 1.0^ of one mA and 
the field produced along the central axis of symmetry was 
determined to within t 10 .% of a nT . 

The data obtained on the sensor's sensitivity for a 
5 nT and 66.4 nT signal, when compared to the theoretical 
curve, revealed that a more uniform field in the vicinity of 
the sensor was needed and that the test should be performed 
in an area remote from metallic structures such as in a 
laboratory. The data obtained was within 6 dB of the theo- 
retical values as can be seen in Figure 12. 

METHOD II. 

The second method, that of using a pair of Helm- 
holtz coils, produced a more uniform field over a larger 
area. Figure 11 illustrates how the coil was designed to 
produce an easily accessible, almost homogeneous magnetic 
field in the central region between the coils. 

Expressions for the magnetic induction vector for 
the axial and radial componets are: 
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54 z^r'^ ^ 
125 R*" 



..] (3-10) 




(3-11) 



as given in (Jefemenko, 1966). 



So, in the central region: 
g ^ Uq3NI 

^ 5/5 R 



(3-12) 



= 0 (3-13) 

to the terms of the order (L/R)^, where L is the linear di- 



mensions of the region. Using the same method to determine 
the current in the loops, as explained above, the data 
shown in Figure 13 indicates that the actual sensitivity 
was about one dB less than the theoretical value for the 
entire frequency band of interest. It was determined that 
the accuracy of the measured sensitivity was well within 
the expected theoretical limits. The major source of error 
in the theoretical computuation was determining accurately 
the area enclosed by the conductor due to the finite width 
of the coil. 
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Figure 9 Magnetic Field Along Axis 
of a Current Carrying Ring 




Figure 10 Hookup for Establishing a Known Field 




Figure 11 Helmholtz Coil Schematic 
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To obtain the corrected power spectral density referenced 
to 1 nT^/Hz, the following procedure (derived in Appendix 
D) was used: 

N dB = The Spectrum Analyzer reading in dB volts. 

(A) = The Hanning windov; bandv/idth correction for 

the particular Spectrum Analyzer used. 

(B) = The transfer Function defined above. 

(C) = The corrected power spectral density where 

0 dB is referenced to InT^/Hz. 

(C) = N dB - (B) - (A). 
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0 dB Re: 1 Volt 
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frequency 




I 
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Fig\)re 13 



C. SYSTEM TRANSFER FUNCTION 



The system can he represented by a simple block diagram 
that converts a varying magnetic field into a voltage that 
can be displayed and whose power spectral density can be 
measured with a Spectrum Analyzer: 



1 nT 



SYSTEM 



X number of volts 



The transfer function is frequency dependent. It was cal- 
culated by separating the system into its components as 
follows: 



1 nT 




Once the sensor sensitivity was determined using the Helm- 
holtz coils it remained to compute the gains and conver- 
sions of the electronic components to determine the overall 
system transfer function. 

The preamplifier was designed to be a low pass, low 
noise amplifier with a very sharp cutoff frequency at 20 Hz. 
This cutoff frequency is necessary to attenuate the 60 Hz 
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fields present at all land sites. As can be observed in 
Figure 14 the gain characteristics are flat to 10 Hz with 
a gain of 66.5 dB. The sharp cutoff at 20 Hz was obtained 
with a final stage active filter network which added slight- 
ly to the gain prior to cutoff. 

The voltage controlled oscillator (VCO) was designed 
to have a center frequency of 1500 Hz with an output voltage 
level of .4 volts (Peak to Peak). The purpose of having a 
low output voltage at 150O Hz is to minimize crosstalk be- 
tween tape recorder channels and to operate well within the 
frequency response of the tape recorder. Since the tape 
recorder and mixer (described in Section III. A) have unity 
gain, the overall gain of the VCO and frequency-to -voltage 
converter (FVC) was measured by amplifying the I500 Hz 
signal to the required level to drive the FVC with a lov/ 
noise, PAR-1 I3 amplifier. The low pass gain characteristics 
of the VCO/FVC combination can be seen in Figure 15* The 
response is similar to that of the preamplifier and is 
characterized by a constant gain of 3^*7 dB from dc to 4 Hz 
and a slight decrease after 10 Hz as required for additional 
60 Hz filtering. 

The system transfer function was calculated by adding 
the gain of the preamplifier to the gain of the VCO/FVC com- 
bination and to the sensor sensitivity at equivalent fre- 
quencies. This can be seen in Figure I6. The data presented 
in this research is referenced to 1 nT rather than 1 volt. 
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Figure 16 



D. SYSTEM NOISE 



The best amplifier performance is realized when the 
amplifier causes the least decrease in the overall signal- 
to-noise ratio. In many instances the thermal noise gener- 
ated by the signal source resistance is the dominant factor 
in determining the input signal -to-noise ratio. Amplifier 
noise performance can be specified by the amount of noise 
the amplifier adds to the amplified source thermal noise. 
Expressed in decibels, and defined as the "noise figure", 
or NF, it is given by; 

NF = 20 log total rms noise voltage at amp output 
gain X (source thermal noise voltage) 

where the Source Thermal Noise = V^KTRAf volts rms, and 
K = Boltzmann's constant, l.Jd x 10”^^ joules/°K 
T = absolute temperature in degree Kelvin 
f = Equivalent Noise Bandwidth in Hz 
R = source resistance in ohms 
The total output noise may be converted to an equivalent 
input noise by dividing by the amplifier gain. Thus the 
noise figure becomes: 

NF = 20 log total rms noise voltage referred to amp input 
source thermal noise voltage (rms) 

The equivalent noise bandwidth (Af) used in determining the 
source thermal noise is simply tr/2 times the high fre- 
quency rolloff. 
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The system's preamplifier has a 3 dB point on the low 
pass gain curve at 25 Hz. This yields an equivalent noise 
bandwidth of 



Af = 25 X it/2 = 39.27 Hz 

and a source thermal noise of 8.69 x 10 “^ rms volts, which 
is -I 6 I dB re: 1 volt. 

A 1201^ resistor was connected across to the preamplifier 
input while the rms voltage was measured on the spectrum 
analyzer. By dividing the total output noise by the amp- 
lifier gain the Noise Figure at one hertz was obtained as 
15 dB (see in Figure 17). 

Once the total rms noise voltage referred to the pre- 
amplifier’s input was determined, the noise introduced by 
the remainder of the system, was determined to be negligible 
in comparison to the preamplifier noise, with the exception 
of four very discrete frequency components introduced by 
the tape recorder. These are centered at 3*^» 6.1, 6.8, 
and 8.1 Hz and have been deleted from the data curves. The 
purpose of the reference signal and mixer is to eliminate 
the tape recorder noise. However, the noise reduction 
is limited to about 10 dB at these four undesired frequen- 
cies . 

The total system noise rms voltage level was measured 
several times at various locations. In particular, it was 
measured each time data was taken in the environment pre- 
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sent (i.e. land or sea). On each occasion the total noise 
was found to he at the level shown in Figure 1? which was 
exclusively front end preamplifier noise. No other sources 
introduced noise above this level. 
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Figure 1 ? 



IV. EXPERIMENTAL RESULTS 



A . UNDERWATER DATA 

Of the six hours of total recorded data, half was ob- 
tained at sea and half at a remote land site near Chew's 
Ridge in California . The underwater data was taken at two 
stations designated as Station #1 (36°-38'N, 121°-52'W) and 
Station #2 (36°-39'N, 121°-53'W). The depth, local time, 
date, and frequency band of the measurements are annotated 
on each figure. The sea floor at both stations has less 
than 5° of slope and consists mainly of smooth sand. The 
sensor was arranged so that the central axis of symmetry 
was parallel to the sea floor at an arbitrary angle refer- 
enced to a North/South direction. 

Figure 18-19 represent one half hour of recorded under- 
water data taken in a moderately heavy sea state ( 10-15 ft. 
swells) on 14 April 80, The remainder of data taken on 14 
April was not usable due to loss of electrical continuity 
between upper and lower spheres. Figures 20-29 represent 
2.5 hours of recorded underwater data taken in calm seas 
on 18 April 80. 

B . LAND DATA 

All three hours of land data was recorded on 1 May 80. 
The location was chosen because of its remoteness and lack 
of external power sources or power lines. The position of 
the sensor remained fixed throughout each data run, and 
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the orientation was as follows: 



VERTICAL 




Figure # 
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70° 


180° 
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NA. 
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90° 


90° 
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90° 


0° 
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60° 


0° 


36 


30° 


180° 


37 


90° 


90° 



The only prominent signals observed at less than 1 Hz were 
on the first run. All subsequent runs revealed preamplifier 
noise to be dominant at less than 1 Hz. 
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C. DISCUSSION AND CORRELATION OF DATA 

Many variables influence the level of geomagnetic acti- 
vity on the sea floor, thus, an unambiguous comparison with 
the sparce previous data is almost impossible. Three at- 
tempts at sea of three hours duration each were made with 
the coil's axis of symmetry directed upward in order to mea- 
sure the vertical field component. The data gathered indi- 
cated that the geomagnetic activity in the vertical plane 
was below the noise level of the system. This agrees with 
theory which predicts that incident plane waves are refracted 
normal to the surface leaving only a component of the field 
in the horizontal plane in the sea. The data indicates that 
a small verticle component in the sea may still be present 
but was observed to be below preamplifier noise. Thus an 
upper limit on the vertical component of magnetic field 
fluctuations can be established as 10“3nT^/Hz at 1 Hz and 
10-6nT2/Hz at 10 Hz. 

It is also difficult to compare the fluctuations of the 
total magnetic field ats obtained from a total field measuring 
system and a system which measures only one component of the 
field. Earlier data obtained with a Cesium magnetometer 
(Chaffee 1979) on the floor of Monterey Bay indicated that 
the 1 Hz fluctuations as measured in a 5 hour period were 
-30 dB referenced to 1 nT^/Hz at a depth of 70 meters. The 
average value obtained with the coil on 18 April at 1 Hz was 
-28 dB referenced to 1 nT^/Hz at a depth of 45 meters. This 
appears to be in very good agreement. 
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The most interesting frequency range observed was above 
1 Hz. The level of geomagnetic activity was observed to be 
quiet between 1-6 Hz but became active at 9*6, 11. 9» 16.5» 
and 19 Hz for the sea measurements. The land measurements 
indicated similar characteristics with predominant fre- 
quencies centered at ?.3> and 24.9 Hz. Polk and 

Fitchen (1962) made measurements of the N-S, E-W horizontal 
magnetic components and showed that the fundamental mode of 
Schumann resonance is not constant in frequency with time. 
They were able to determine the spectral frequency range as 
7 Hz to 11 Hz. The spectral frequencies at 7*3 Hz and 9*6 
Hz as measured on land and in the sea respectively, corre- 
late well with the fundamental Schumann resonance frequency. 

The spectral frequencies of 10. 7 Hz and 11.9 Hz as 
measured on land and in the sea respectively have the typical 
time development for Schumann resonances as reported by 
Polk and Fitchen. The second resonant mode should occur at 
14 Hz as observed by Balser and Wagner (I960). Since the 
data in this research does not correlate well in frequency 
as a second resonant mode, but does indicate typical wave- 
shapes for Schumann resonance, furthur investigation is 
recommended. 

The 7*3 and 10.7 Hz spectral lines appeared to have 
different polarizations as different coil orientations pro- 
duced varying levels of either frequency. This observation 
could not be established on the sea floor, as no control 
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over the angle with respect to North/South was available, 
but was quite evident in the land measurements. The 24. 9 
Hz spectral line was observed to increase toward the after- 
noon and was apparently independent of orientation. Fur- 
thur investigation of these frequencies is required to de- 
termine polarization as the methods utilized for determining 
orientation were inaccurate (t 5°) but the data was still 
reproducable and quite convincing. 

The sea floor data taken on 14 April was higher in 
magnitude than the data taken on the 18 April. In part this 
can be attributed to a higher A-index (9) on 14 April com- 
pared to an A-index (4) on 18 April. Internal and surface 
wave motion was higher on 14 April as evidenced by separa- 
tion of the cable. No wave height correlation could be 
made, since a wave height measuring device was not avail- 
able in the vicinity of the data gathering region. 
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V. equipment/system improvements and recommendations 



The system utilized for this project was well suited 
to a shallow sea water environment (200 meters). The 
system and equipment have been demonstrated to operate as 
conceived. To deploy this system in a deep ocean environ- 
ment will require several modifications. 

The following discussion on system improvements is 
presented as recommendations to make the system more reli- 
able and usable in the ocean environment. The tape recorder 
is probably the single, major source of equipment troubles 
in the system. Although the record heads and capstan drive 
were meticulously cleaned and demagnetized, the recorded 
information still appeared noisy and distorted. The mixer 
circuit was able to attenuate and filter some of these 
distortions, but not to the degree that totally eliminated 
them. Initial investigation revealed that the longer the 
tapes, the more prominent the noise. Shorter tapes were 
utilized with some success. Further investigation revealed 
that the single capstan drive motor and associated circuitry 
were of inferior design and failed to adequately compensate 
for the changing tape reel speeds. As discussed in a pre- 
vious section, the distortions were finally determined to 
be present at certain discrete frequencies. 

Another obvious requirement for a more reliable re- 
corder is the reduction of recording time brought about 
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by the shorter tapes. Even though digitizing the data 
would increase data gathering time, the present recording 
system's operation in this role is highly questionable. 
Therefore, it is recommended that this tape recorder be re- 
placed by a more reliable instrument. 

Coils used as sensors require less power to operate 
than the more complicated total field measuring magneto- 
meters. However, coil type sensors have lower sensitivities 
(sensitivity is proportional to frequency) and are unable 
to measure the dc component of the earth's main field. The 
present coil size is restricted by the largest glass sphere, 
as discussed in Section III. The sensitivity of the sensor 
can be increased by increasing the area and/or the number 
of turns of the coil. However, this is not possible in the 
present system. Thus the sensor is limited in its capa- 
bilities to detect geomagnetic fluctuations because of its 
lack of sensitivity at low frequencies. 

In addition, higher frequencies are attenuated more 
than lower frequencies in sea water (see Section II). 

Hence, the use of a coil in a deep ocean environment is 
questionable but not completely ruled out. Modifications 
to the preamplifier (to be discussed below) would increase 
the sensitivity of the system and, possibly make it quite 
feasible for deep ocean use. However, it is recommended 
that continued use of the present system be restricted to 
shallow water deployment or the system be redesigned for 
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land measurements. That is, a larger sensor be designed and 
constructed for increased sensitivity and used to gather 
geomagnetic fluctuations on land. If three such sensors 
were produced, then all components of the earth's magnetic 
field fluctuations could be measured. In addition, the 
polarization of the horizontal component could be further 
investigated by connecting two such coils and determining 
the coherence between the two measurements. 

State of the art devices limit the design of a pre- 
amplifier which operates at low frequencies. In the present 

system, the noise figure of the preamplifier is 15 dB at one 
hertz. More recently an improved preamplifier has been 

designed and built by Dr. Alan Phillips whose noise figure 

is about 10 dB better. This improvement will greatly 

enhance the sensitivity of the system and make its use in 

a deep ocean environment more feasible. 

In the present system, no means exist for determining 
the horizontal orientation of the sensor axis when deployed 
in the ocean. Observed data indicates that the horizontal 
component of geomagnetic field variations may be polarized. 
(See Section IV on discussion of data). It is recommended 
that an aiming device be designed and constructed so the 
exact orientation of the sensor is known when deployed in 
the ocean. 

An improvement in the splicing of cable and connectors 
is recommended. Presently the Research Department at 
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the Naval Postgraduate School splices the connector pigtails 
to the cable and sphere penetrators. This method was 
satisfactory for shallow water employment, however, a more 
reliable method is required for deep ocean usage. It is 
recommended that cables, connector pigtails, and penetrators 
be commercially spliced for greater reliability. 

Presently the only means of data correlation are with 
past data collected at the Naval Postgraduate School and 
the indices (see Section II) reported by Fredericksburg, 
Virginia. Since surface and internal ocean waves are known 
to be a source of field fluctuations, it is recommended that 
a wave gauge (various types are available) be purchased and 
utilized as another means of correlating geomagnetic 
variations . 

It is recommended that the various electronic component 
chassis's be integrated into one or two units and a more 
practical method of connecting power be provided to en- 
hance the system’s compactness and increase its reliability. 
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APPENDIX A 



EQUIPMENT SCHEMATICS 

Numerous preamplifiers were tested and utilized during 
the extent of this project. Figure 38 was the most recent 
"State of the Art" amplifier which provided the quietest 
operation, i.e., greatest signal -to -noise ratio. 

Figure 39 is the schematic for the mixer circuit used 
at the output of the analog cassette recorder for data 
analysis. Figure 40 is the 2 KHz oscillator reference 
signal genera.tor used in the analog recording system to 
reduce tape "wow and flutter" . Figure 4l is the VCO unit 
used in this experiment. 
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Figure 4l 



APPENDIX B 



EQUIPMENT USAGE 

Equipment utilized to measure, record, and analyze 
magnetic field fluctuations requires a certain amount of 
maintenance and handling. Below is a suggested list of 
hints which might help prevent equipment damage and pro- 
long its usefulness o 

A. Benthos glass spheres - much care should be exer- 
cised when handling the hemispherical shells because they 
are highly susceptible to chipping around the groimd glass 
edges when exposed o The shell halves are serialized and 
appropriately marked for proper mating. When exposed to a 
deep water environment, the spheres are subjected to lower 
temperatures. Retreaval from the ocean depths results in 
the condensation of water. A heat gun may be necessary to 
dry the surfaces and return the spheres to ambient temper- 
ature. Temperature differences may also cause the glass 
housings to become vacuum locked. Evenly applied heat 
about the sphere surface will assist in releasing the 
vacuum. 

The strip sealant used for sealing the shells 
should be applied on a meticulously cleaned surface at or 
near room temperature. Failure to properly seal the 
spheres will result in leakage and possibly, equipment 
damage. Toulene and alcohol are the recommended glass 
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cleaning solutions. Caution; Toulene is both toxic and 



poisonous. It should be applied with rubber gloves in a 
well ventilated area. 

B. Electronic Equipment - of particular concern here 
is proper application of power supply voltages. Improperly 
applied voltages may result in integrated circuit and/or 
printed circuit board damage. The suggested application of 
dc power is to first, connect the non-common lead (s), 
being careful not to short the common lead against the 
equipment chassis, and finally, connect the common lead. 
Circuit power should not be applied unless there is a load 
on the input . 

Co Portable Cassette Tape Recorder - The tape recorder 
heads and capstan drive should be cleaned and demagnetized 
after each recording. Failure to do this may result in 
noise spikes that could ruin a hard days work.. It is also 
recommended that all tapes be "fast forwarded and rewound" 
before being utilized for data taking. High quality 60 
minute tapes are recommended, since longer tapes appear to 
impose greater drag on this type recorder. Before sealing 
the instrumentation sphere, ensure the record button is 
depressed and the record levels are properly adjusted. 

D. Voltage controlled oscillator (VCO) - the center 
frequency of the VCO should be adjusted to I500 Hz prior to 
sealing the instrumentation sphere. This can be accomplished 
by adjusting the center frequency poteniometer with the input 



shorted. 
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E. Mixer - an amplifier stage is provided for the 
sensor and reference oscillator (clock) inputs to the 
mixer. With inputs from the clock and sensor, adjust the 
gain poteniometers for a 6-10 V peak-to-peak output dif- 
ference frequency. A 6-10 V peak-to-peak output is 
necessary to drive the frequency-to-voltage converter. 

Verify that the output difference frequency is what it 
should be. Note that this procedure may have to be per- 
formed prior to each data analysis, since record levels 
may vary from one recording to the next . 

F. Gelyte batteries - the batteries should be charged 
prior to taking data. Each battery is labeled and a record 
of charge is maintained in the laboratory. Presently there 
are 2 ampere-hour ratings being utilized in the system. 

The recommended charge procedure is: 

18 ampere-hour ; charge at 1 ampere for approxi- 
mately 20 hours for full charge. Do not overcharge for more 
than 24 hours in excess of the ampere-hour rating. 

2.6 ampere -hour ; charge at 0.2 amp for approxi- 
mately 15 hours for full charge. Do not overcharge for more 
than 24 hours in excess of the ampere-hour rating. 

G. Connectors/ cable ; before and after each use in a 
corrosive environment, all connectors should be scrupulously 
cleaned to prevent deterioration. Alcohol and cotton-tipped 
applicators should be used. Further, all male connectors 
should be lightly lubricated with silicon grease to ensure 
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an adequate seal and prevent water intrusion. The locking 
sleeves on the connectors are there to prevent stress and 
are not part of the pressure system. The connectors are 
adequately marked and care should be exercised in mating 
them. The cable and connector pigtails should not be sub- 
jected to undue strains and sharp bends. Intermittent 
and/or complete signal loss may result as well as cable 
separation. 
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APPENDIX C 



Geomagnetic Indices and Forecasts can be a valuable aid 
in data correlation and interpretation. In this Appendix 
we include an explanation of the type of geomagnetic in- 
formation which is available on a weekly basis from NOOA 
(Boulder, Colorado) and which we have found to be very 
helpful in our research. 
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Basking in its ixu'^tk, feeding on its energy, 
plagued by its amptiens, man ponders the 
mysteries of the sun. 



bVAA-SESC photograph of Galen McFOiiyen atop 
Boulder *s Flatiron fountains by Joe Sutorik, 
December, 1968. 
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DESCRIPTIVE TEXT 



Concents of 

Preliminary Report and Forecast of 
Solar-Geophysical Activity 

1 May 1980 



The Preliminary Report and Forecast of Solar-Geophysical Activity 
(PRF) is composed each Tuesday and issued the next day by the joint 
National Oceanic and Atmospheric Administration, Space Environment 
Laboratory (N0AA)/Air Force Air Weather Service (AWS) , Space 
Environment Services Center (herein denoted as SESC) , Boulder, 

Colorado 80303. It is a continuation of the series of weekly 
retorts which began in 1951 and were issued jointly by the High 
Altitude Observatory and NOAA and its predecessors. The current 
series began with serial number PRF 001 on 09 September 1975. 

This is a preliminary reoqrt based on data available at the time 
of ^oblication and, therefore, cannot reliably be cited for reference 
purposeT ^ It is intended for rapid distribucion and ox primary use 
to commercial and real time operations. Many of the data contained 
here are subject to later revision or refinement. There are no 
copyright or other restrictions in photocopying this report. A.rchival 
quality data, suitable for more extensive studies, are published in 
the monthly publication ’’Solar Geophysical Data” (SGD) available 
from the National Geophysical and Solar-Terrestrial Data Center 
(World Data Center A), Environmental Data and Information Service, 

D6 NOAA, Boulder, CO 80303. Ac a still later date, definitive 
publication for much of the data is made in the lAU Quarterly 
Bulletin of Solar Activity. Special reports are also issued from 
time to time by the World Data Center for Solar-Terrestrial Physics. 

This report is made possible through the combined efforts of 
several groups providing data on a fast time basis. A complete 
list of data sources is given in Note 4 at the end of this descrip- 
tive text. Comments concerning the content of or possible improve- 
ments in Che Report are invited. 

Frontispiece 

The frontispiece contains the title of the publication, publication 
number and date, and a photograph symbolizing man's dependence upon 
the radiative output of the sun. The photograph is a frame from 
the time-lapse 35 mm patrol movie taken with the hydrogen-alpha 
patrol telescope of the Space Environment Services Center in Boulder, 
Colorado. The photograph was taken December 9, 1968 by Joe Sutorik, 
chief solar observer, in collaboration with student observer and 
mountaineer Galen McFadyen. The two observers carefully charted 
the positions of the setting sun along the profile of the rugged 
Flatiron Mountains three miles west of the observatory and antic- 
ipated Che date when the sun would set behind an accessible rock 



92 



outcropping. On the selected date, Mr. HcFadyen scaled the mountain 
and positioned himself so that his full figure would appear silhou- 
etted against the sun. 

First pag;e 



The first page summarizes the highlights of solar-geophysical 
activity during the past seven days and contains forecasts of 
expected activity during the next 27-day interval. 

In the activity summary and forecasts, solar active regions are 
identified by the SESC-assigned region number, latitude, Carrington 
(heliographic) longitude, and modified Zurich region classification/ 
sunspot area in millionths of the solar hemisphere on the date of 
maximum sunspot area. 



Significant solar activity is discussed in terms of the character- 
istics of the region of origin, x-ray classification of flares (C, 

M, or X), optical importance of flares (SN, IB, etc.), radio emission, 
energetic particle emission, and geophysical effects. The character- 
istics of near-Earth proton events as detected by both satellites 
and ground-based sensors are discussed. Geomagnetic storms and 
disturbances are described. Observations of visual aurora are 
described when available. Whenever feasible, solar and geomagnetic 
activity will be summarized and forecast within the framework of 
the following routinely used terms: 



Solar activity: 



Very low 



Low 



Moderate 



High 



usually only quiet regions on the solar disk 
and no more than five of these; fewer than 
ten C-class subflares without important 
centimetric radio bursts or sudden ionospheric 
disturbances during a 24-hour period. 

usually more than five but less than ten 
quiet regions on the solar disk; only C-class 
subflares without important centimetric radio 
bursts or sudden ionospheric disturbances 
during a 24-hour period. 

eruptive regions on the solar disk; fewer 
than five M-class x-ray events with important 
centimetric radio bursts or sudden ionospheric 
disturbances during a 24-hour period. 

active regions on the solar disk; several M- 
class x-ray events with important centimetric 
radio bursts and strong sudden ionospheric 
disturbances and/or one or two importance 2 
chromospheric flares or X-class x-ray events 
during a 24-hour period. 



B. 
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Very high - one or more regions on the solar disk capable 
of producing proton events; one or more 
chromospheric flares of importance 2 or 
greater with outstanding centimetric radio 
bursts (500 flux units above background, or 
greater), X-class x-ray events, and major 
sudden ionospheric disturbances observed or 
expected within a 24-hour period. 

The x-ray classification of solar flares was initiated by the Space 
Environment Services Center on 01 January 1969. This classification 
attempts to classify solar activity by its energy output, and 
offers at least two distinct advantages compared to the standard 
optical classifications. It gives a better measure of the geophys- 
ical significance of a solar event, and it provides an objective 
means of classifying almost all geophysically significant activity 
regardless of its location on the solar disk or at (or beyond) the 
solar limb. The definitions of class C. M, and X events are based 
upon the peak x-ray flux in the 1 to 8 Angstrom wavelength range, 
as follows: 

Class Peak l-SA flux 

in watt m’^) 

C $ < 10"^ 

M 10‘^< -t < 10“^ 

X * > 10"^ 

For descriptive purposes, a one digit number from 1 to 9 may be 
appended to these letter designations. The letter then acts as a 
multiplier. For example, a C3 burst indicates^an x-ray burst with 
peak flux in the neighborhood of 3x10"^ watt ra Similarly^ an X5 
event would indicate an x-ray burst with peak flux near 5x10 ^ watt 
m** . This expanded classification is generally used only for 
description of events; forecasts are usually issued in terras of the 
broad C, M, and X categories. CO^designation may be used for 
events smaller than lxl0”6 watt m“2. 

Geomagnetic activity: 

In the definitions of the following adjectives used to describe 
geomagnetic activity, A refers to the 24-hour A-index observed at a 
mid-latitude observatory such as Fredericksburg, VA; not the planetary 
A-index, Ap . The K-indices are, likewise, mid-latitude values. 



Quiet 


= A 


<_ 7, usually no K-indices > 2. 


Unsettled 


= 7 


< A <15, usually no K-indices > 3. 


Active 


= 15 


< A < 30, a few K-indices of 4. 


Minor storm 


= 30 


< A <50, K-indices mostly 4 and 5 


Major storm 


= A 


> 50, some K-indices 6, or greater. 



-2 -1 
(o in erg cm sec ) 

$ < 10 "^ 

lo'^l « < 10'^ 

$ > 10 '^ 



c. 
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Sudden commencements are indi.cated by beginning times given to the 
nearest minute. Gradual commencements are indicated by beginning 
times to the nearest hour. 

Immediately following the summary of activity and the forecasts is 
a computer-generated display of classical solar-geophysical indices 
for the dates indicated. The column headings are defined as 
follows : 

OTT 2800 MHz: The Ottawa NRC 10 centimeter (2800MHz) radio flux 
observed at approximately 1700Z on that date. 

SSP NBR: The official SESC sunspot number reported for that date in 
the Joint USAF/NOAA Solar Region and Activity Summary issued at 
0200Z the next day by SESC. Sunspot reports are received daily 
from as many as half a dozen observatories. The SESC staff 
reviews these reports and forms a composite picture of each 
individual region, taking into account such factors as the time 
of observation and the quality of seeing. The composite informa- 
tion for all regions is then incorporated in the 0200Z message. 
The sunspot number is computed according to the formula R=10G+N 
where G is the number of sunspot groups (regions) and M is the 
total number of individual spots comprising all groups. A 
sunspot number of zero indicates there were no visible sunspots 
on that date, v;hile a blank indicates that no observations 
V7ere possible. 

Fredericksburg and Anchorage A- and K-indices: The daily 24-hour A- 
index and 3-hourly K-indices from the Fredericksburg VA (mid- 
latitude) and Anchorage AK (high- lat itude) stations monitoring 
the Earth’s magnetic field. K-indices range from 0 (very 
quiet) to 9 (extremely disturbed). A-indices range from 0 
(very quiet) to 400 (extremely disturbed) with an A-index of 
30 or greater, indicating geomagnetic storm conditions. 

Additional pages 

Additional pages are included in the Preliminary Report as solar 
activity requires, and are comprised of computer-generated displays 
which are described below in their normal sequence. The sequence 
and page numbering may vary occasionally depending upon the level 
of activity, reports of special events, or the inclusion of data 
for past days which could not be published at the usual time. 

Alerts Issued This Week 

A summary of all events which exceeded the SESC alert thresholds 
during the previous week. Listed are the dates, times and the 
activity that resulted in the issuance of the alert. This summary 
will enable alert customers to have a complete listing of alerts 
and may adequately fill the needs of seme customers currently 
contacted by telephone. 



D. 
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Weekly Energetic Event Summary 

A suTTtinarv of all_energetic (x-ray intensity greater chan or equal 
Co lxlO~^ watt m”2) events occurring during the seven days preced- 
ing the Report composition date. , Reported are: 

Month and Day: Date event occurred. 



Begin, Maximum and End: Begin, maximum and end times of 

the x-ray burst. The begin time is often well defined 
due to the impulsive beginning of most flare associated 
x-ray events. In the past, x-ray events were ended 
when the flux levels decayed to a point halfv/ay between 
the maximum level and the pre-flare background level. 

This definition occasionally resulted in unterminated 
events. For operational purposes the end time of an x- 
ray burst is now determined by both the peak flux and the 
decay rate. An X-class event is ended when flux levels 
decay to less than M5 levels and remain stable for a 
minimum of four minutes or when flux levels decay to less 
than the Ml level, regardless of stability. All bursts 
smaller than X-class are now ended when flux levels decay 
to the C5 level or when flux levels decay to a point 
halfway between the maximum level and the pre- flare 
background level and remain stable for a minimum of four 
minutes. The typical x-ray signature dictates that the 
end time of a burst will be less well defined than the 
begin time. 



X-ray Class: As defined earlier in this descriptive text. 



OPT; Optical information (if available) as follows: 

Type of Activity: FLA - Solar flare 

EPL - Eruptive prominence on the 
solar limb 

BSL - Bright surge on the limb 
SPY - Spray (including ejecta) 

LPS - Loops 

DSF - Disappearing solar filament 
UNK - Unknown activity 
Importance; S = Subflare 

area < 2.0 square degrees 

1 = Importance 1 

2.0 _< area £ 5 . 1 square degrees 

2 = Importance 2 

5.2 _< area £ 12.4 square degrees 

3 = Importance 3 

12.5 ^ area _< 24.7 square degrees 



4 = Importance 4 

24.8 <_ area in square degrees 



E. 
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Brightness qualifiers F, N or B are generally appended to 
the importance character to indicate Faint, Normal or 
Brilliant, respectively. 

Location: Flare location in heliographic coordinates. 

Region: SESC-assigned region n\mber. 

Radio: Flux unit peak value above pre-burst background 

of associated radio bursts at frequencies of 245, 

415, 606, 1415, 2695, 8800 and 15400 MHz. 

1 Flux Unit = 10”^^ watt m"^ Hz”^. 

Sweep: Sweep frequency radio event associated with the energetic 

event, as follows (generally only II, IV and CTM are reported: 



Type; 


I 


Noise storm 




II - 


Slov7 drift burst 




III - 


Fast drift burst 




IV - 


Broad-band smooth continuum burst 




V 


Brief continuum burst generally associated 
with Type III bursts 




CTM - 


Long-lived continuum associated with storms 
of Type IV bursts 




STM - 


Storm of Type III bursts 




MWB - 


Short microwave continuum 




SDF - 


Decimeter slow drift burst 


Importance : 


Relative importance on a scale of 1 to 3. 

An importance 1 sweep event would be considered 
small or minor. An importance 3 event is a 
major or significant event. 



SID: Sudden Ionospheric Disturbance information, 
as follows: 

Type: SWF - Short-wave fadeout 

Frequency: Highest frequency Sl‘TF was detected. 

Importance: Subjective measurement of the relative 

importance of the event on a scale of 
1 to 3 (much the same as already described 
for sweep frequency radio events) . 

Flare List 



A listing of solar flares (optical importance and/or x-ray class) 
reported during the previous seven days. The list has been edited 
on a daily basis by the duty solar forecaster to delete multiple 
reports of the same event. Criteria used in determining which 
report, or combination of reports, will be used include: Quality of 
the observation (observing conditions, was flare observed in its 
entirety?) , type of report (real time or film review) , and the 
consensus of reporting observatories. The list includes the following 
information: 



R 
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Month and Day: Date event occurred. 

Begin, Maximum and End: Begin, maximum and end times 

of the optical flare. If the start of the flare was not 
observed, a "B" is appended to the begin time indicating 
the flare started "Before” that time. If the maximum was 
not observed, the space is left blank. If the end of the 
flare was not observed, an "A” is appended to the end 
time indicating the flare ended "After” the end time. 

Location: Flare location in heliographic coordinates. 

IMP: Importance of the optical flare, as defined earlier 
in this descriptive text. 

REG: SESC-assigned region number or left blank for unassociated 
x-ray events and for flares which occasionally occur in 
undefined regions. 

Class: Denotes x-ray class of the event according to 

the C, M and X classification described earlier in 
this descriptive text. 

Weekly Region Summary 

A summary of all regions visible on the solar disk during the 
previous seven days. Each date's assessment of a region is based 
upon observations made during the 24-hour UT day and the location 
is corrected to 2400Z on the date indicated. 

This summary information is compiled from the Joint AFGWC/SESC 
Solar Region and Activity Summary which is a composite picture of 
each region from all sunspot reports received at Boulder during the 
UT Day. The Weekly Region Summary contains the following information: 

MO and DA: Month and day 

Region Location: Location of the "midpoint of the latitudinal 

and longitudinal extremeties" of the white-light sunspot 
group in heliographic coordinates at 2400Z on the date 
reported. The position of the center of H-alpha plage 
associated with the region is given if sunspots were net 
visible. 

Sunspot: Sunspot information is as follows: 

Area: Total area of the region's sunspots in millionths 

of the solar hemisphere. 

Z: Modified Zurich (McIntosh) three letter (Zpc) 

classification as follows: 

Z = Modified Zurich Brunner classes, A through H 
p = Penumbra type of largest spot in group: 

X - no penumbra, r - rudimentary, s - small 
symmetric, a - small asymmetric, h - large 
symmetric, and k - large asyirnnetric 
c = Relative sunspot distribution or compactness of 

group: X - single spot, o - open, i - intermediate, 
and c - compact 



G. 
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NR: Total number of individual sunspots in the group or 

region . 

MAG: Magnetic classification, as shovm below: 

A = Alpha (single polarity spot) 

B = Beta (bipolar spot configuration) 

G = Gamma (atypical mixture of polarities) 

BG = Beta-Gamma (mixture of polarities in a dominantly 
bipolar configuration) 

BD = Beta with Delta configuration 
BGD = Beta Gamma v;ith a Delta configuration 

D = Delta (opposite polarity umbrae within single 
penumbra) 

H-Alpha Plage: Indicates no sunspots visible but 

region still visible in K-alpha spectro- 
heliograms . 

Calcium: Calcium spectroheliogram data from Manila, Philippines, 

as follows: 

Area: Total region area of calcium plage in millionths 

of the solar hemisphere. 

BRT: Calcium plage brightness on a subjective scale of 
1.0 (low) to 5.0 (high). 

Blank space indicates no calcium plage report available 
for that region on that date. 

Flares: Tabulation of the number of C, M and X-Class x-ray 

bursts, number of subflares and number of importance 1, 2 
and 3 flares observed on that date. 

At the bottom, the average Carrington longitude of the region and the 
x-ray burst and flare totals are listed. Finally, a comment on the 
fate of the region is included. This indicates whether the region 
crossed over the solar limb, died on the solar disk, or was still 
an active region on the disk at the end of the reporting period. 

Boulder Geomagnetic Substorm Log 

This is a tabulation of substorm occurrences as observed by ground 
magnetometers in the North American International Magnetospheric 
Studies (IMS) network. This network is comprised of twenty-six 
stations; extending from north of the auroral oval to the equatorial 
zone. A substorm is a localized geomagnetic disturbance which 
usually occurs near local midnight and is restricted in longitude. 
However, the current systems developed during a substorm may affect 
ground magnetometers in the entire nighttime sector. Additionally, 
individual substorms may occur at local times, away from midnight, 
and may be as large as 24 hours (global) in longitudinal extent. 

By noting the time, location and scale of a substorm; one may 
estimate the effect of a substorm at a specific location. Among 
the many substorm effects are ionospheric effects (which influence 
radio communications) and telluric effects (which may disturb long 
distance electric power and communications systems) . The familiar 
aurora is a visible manifestation of the geomagnetic substorm. 
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Currently, the Log provides the date, onset time (UT) and direction 
(from Boulder) of each substorm. The direction is listed as '’east", 

’’west” or ’’centered” (at Boulder's geomagnetic longitude). The 
comment section further describes the geomagnetic field for a 
particular day. Occurrences such as sudden commencements (sc), 
sudden impulses, magnetic storm periods, multiple substoirms and 
substorm (relative) strengths are included. 

Recent Monthly Solar Indices (Preliminary) 

This is a listing of monthly preliminary solar indices for the 
past year (updated v;ithin one week after the end of each month) . 

The daily sunspot numbers appearing on the first page of the 
Report, V7hich are a consensus of observations from as many as half 
a dozen observations, are simply averaged and is reported as "SESC 
Observed” . 

Official Zurich sunspot numbers are reported as ’’Zurich” as they 
become available. The ratio of Zurich to SESC sunspot numbers are 
listed under ’’Ratio Zurich/SESC” ; these can be used to determine 
an SESC ”k” (correction) factor. In order to provide the best 
early estimate of the Zurich sunspot number, e.g. near the end of 
the month when the number is not yet available, an SESC corrected 
number will appear in parentheses in the ’’Zurich” column. 

The monthly mean value of 10.7 centimeter radio flux (1 flux unit 
= 10~22 watt m*2 Hz*'^) is listed for comparison with monthly mean 
sunspot numbers for comparable periods. These radio flux values 
have been adjusted to 1 A.U., except where values appear in paren- 
theses. Parenthetical values are the average of the daily values 
transmitted to SESC in real time. Adjusted values are those which 
appear in Solar Geophysical Data (SGD) . 

Special pages 

To some extent the Preliminary Report is regarded as an experimental 
publication by the NOAA/SEL Space Environment Services Center. The 
Center publishes new data, often on a trial basis; within the 
following limitations: 

1. There must be reason to believe that new data will 

be of interest to a broad range of Pveport subscribers. 

2. Inclusion of nev; data must not have a substantial 
impact on the costs of producing the Report (as 
would, for example, publication of photographs). 

Some items added to the Report are of a ’’one- time-only” nature. 

Among these are special drawings, descriptions or reports of unusual 
solar-geophysical features or events. Other data forms, such as 
the tabulation of monthly mean sunspot numbers, are candidates for 
permanent inclusion in the Report. Because a full discussion of 
special items is not possible in this descriptive text; explanations 
will accompany new data forms. Subscribers are encouraged to 
comment on the utility of various new data provided in the Report. 
Suggestions regarding new items to be placed in the Report will be 
carefully considered. 

I. 
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SMS/GOES X-ray and Particle Plots 

The Synchronous Meteorological Satellites (SMS-1 launched in June 
1974 and SMS-2 launched in February 1975) and the Geostationary/ 
Operational Environir.ental Satellites (GOES-1 launched in October 
1975, GOES-2 launched in June 1977 and GOES-3 launched in June 
1978) carry 0.5 - 4A and 1 - SA solar x-ray detectors and seven 
proton sensors covering the range 0.8 to 500 MeV. These sensors, 
as well as a number of other sensors, are currently monitored 
continuously (except during eclipse periods) in real time by the 
NOAA/SEL Space Environment Services Center in Boulder. GOES-2 and 
GOES-3, in geostationary orbit at 110® West and 135® West, respec- 
tively, are now in use; and SMS-1 and -2 and GOES-1 are in standby 
mode . 

Preliminary graphs of SMS/GOES x-ray and proton data are printed in 
the Report as solar activity warrants, with noteworthy x-ray events 
being generally more frequent than solar protons events. 

The top line of the display identifies the satellite, the sensors, 
and the traces, as follows: 

Satellite: Either GOES-1 or GOES-2 

Sensors: Either XRA (x-ray) or EPS (energetic particles) 

Traces: XRA (line) 1 - SA 

(points) 0.5 - AA 

EPS (line) Flux > 10 MeV 

The second line indicates the date and time of the start of the 
plot. The number appearing in parentheses, either (1) or (5) 
indicates whether the plots are based on 1-mlnute or 5-minute 
averages of the data. 

The ordinate units of the XP^A plots are in watt m and the corre- 
sponding C, M and X x-ray classification threshold is labelled on 
the right margin of the plot. The ordinate units of the EPS plots 
are in P cm 2 sec i ster The abscissa is a time scale in minutes 

or hours, indicating the time elapsed since the beginning of the 
plot at the specified start time. 

Bad data have been deleted wherever possible. 

Last page 



The last page contains the current preliminary H-alpha synoptic 
chart. These maps depict the large-scale distribution of solar 
magnetic fields as inferred from the positions of filaments (cross- 
hatched) , plage corridors (solid lines through stippled areas) , 
filament channels (solid lines outside active regions) and sunspots 
(large dots, dra\/n very schematically). The stippled areas represent 
the H-alpha plages. Dashed lines are interpolations and e.stimates 
required to obtain consistency with polarities and patterns observed 
on previous solar rotations. The numbers identify active regions. 

The ’’plus" and "minus” signs indicate the polarity of the solar 
magnetic field. The heliographic latitude is given along the left- 
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and right-hand side of the nap, the Carrington longitude is printed 
along the botton, and the date of the central meridian passage of a 
given longitude appears at the top. Some portions of the naps may 
be inaccurate; however, every effort is made to correct and refine 
the maps as more detailed study and comparison of tine sequences of 
H-alpha photographs and solar magnetoerams permit. Edited charts 
for comnlete Carrington rotations appear later in Solar Geophysical 
Data (SGD) . 



Notes : 

1. All times used within the Preliminary Report are Universal 
Times or Greenwich Mean Times (as indicated by UT, GMT or 
Z). 

2. The highlights on the first page and a] 1 data listings 
cover the seven-day period Monday through Sunday. The 
27-day forecast period begins on the following Wednesday. 

3. All data listings are prepared by computer. Occasionally, 
inaccurate data will be reported and included. SESC 
appreciates notification of any questionable data. 

4. The following are observatories whose reports are included 
in the SESC data base and utilized for preparing the 
Preliminary Report and Forecast of Solar-Geophysical 
Activity: 
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STATION 

Anchorage, Alaska, USA 
Athens, ‘ Greece 
Boulder, Colorado, USA 
College, Alaska, USA 
Culgoora, Australia 
Ft. Yukon, Alaska, USA 
Fredericksburg, Virginia, USA 
GOES-2 (USA satellite) 

GOES-3 (USA satellite 

Hiraiso, Japan 

Kolloman AFB, New Mexico, USA 

Irkutsk, USSR 

Kitt Peak, Arizona, USA 

Learmonth, Australia 

La Posta, California, USA 

Manila , Philippines 

McMath Obs., Pontiac, Michigan, USA 

Moscow, USSR 

MSFC (Huntsville), Alabama, USA 

Ottawa, Canada 

Palehua, Hawaii, USA 

Poitiers, France 

Ramey, Puerto P%.ico 

Sacramento Peak, New Mexico, USA 

Sagamore Hill, Massachusetts, USA 

Sheep Mountain, Alaska, USA 

Sydney, Australia 

Thule, Greenland 

Tokyo, Japan 

Vos tok, Antarctica 



VJMO ID 


URSl ID 


LOCATION 


70273 


25601 


N61 


W150 


16716 


32401 


N38 


E 23 


72469 


20401 


N40 


W105 




25602 


N64 


W147 




85303 


S30 


E149 


70194 




N67 


W145 


72405 


18403 


N38 W 77 
W109 

geostationary 

W135 

geostationary 




44401 


N36 


E141 


72269 


21305 


N33 


W106 


30710 


40501 


N52 


E104 




21304 


N32 


Will 


94302 


81202 


S22 


E114 


72287 


22301 


N34 


W113 


98429 


42101 


N14 


E121 


72537 


13404 


N42 


W 83 


27612 


34502 


N55 


E 37 




19301 


N35 


W 87 


72628 


18406 


N45 


W 76 


91178 


26204 


M21 


W158 


07335 


30511 


N46 


E 0 


78514 


17201 


N18 


W 67 




21301 


N32 


W105 


72509 


17401 


N42 


W 70 


70269 




N62 


W148 


94768 


85304 


S34 


E151 


04202 


17801 


N76 


W 68 


47662 


44406 


N36 


E140 


89606 


81801 


S78 


E105 
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SmQ^AKY OF SERVICES FROM THE NCA/. SPACE EHVIRC^u•IE^;T SERVICES CENTER 



The weekly publication "Prelin:inary Report and Forecast of Solar 
and Geophy s icaF'TvctTvTty . " A summary of activity indices and 
.event data for the preceding v;eek and a prediction for the next 
27 days. 

A tape recorded message of solar and geophysical activity and 
indices for the most recent 24 hours and for the next 24 hours. 
Telephone numbers are: Coirmiercial : (303) 499-8129 or FTS : 

323-3235. 

A 40-second message with similar information on at 18 

mirmtcs past each hour. 

Primary Report of Solar and Geophysical Acti/itv issued each 
day at 2200UT. Contains activity summary, solar and geomagnetic 
activity indices for previous 24 hours and predictions for the 
next 72 hours. Available by teletype v.^ith recipient expected 
to pay cost of transmission. 

HF Radio Propagation Report issued each day at 0600UT (with 
secondary messages at OUOOUT, 1200UT, ISOOUT) . Distribution is 
also by teletype with recipient expected to pay costs. 

Direct inquiries to the duty forecaster (0700 to midnight, seven 
days per v;eek, Mountain Time) . Call (303) 499-1000 extension 
3171 or FTS 323-3171 (We do not accept collect calls) . 

Event notificatio n by phone or Teletype. Persons or organizatic 
requiring jimmediate notification of the prediction or occurrence 
of various solar geophysical phenomena are contacted v/hen their 
alert thresholds are met. 

Direct access to the SELDADS (Space Environment Laboratory Data 
Acquisition and Display System) . Using standard computer termin 
customers may access the SELDADS and obtain printouts of solar, 
interplanetary and geomagnetic variations data. 

Descriptive leaflets of each of these are available. For more 
information, please wnrite to: 

Space Environment Services Center 
NOAA R43 
325 Broadway 

Boulder, Colorado 80303, USA 

The Space Environment Services Center is a real time service 
organization for serving those people needing real time or 
near-real-time operational support. If you need archival data 
for research studies, you should write to: 

National Geophysical and Solar-Terrestrial Data Center 
Environm»ent Data Services 
D6 NO.AA 
325 Broadway 

Boulder, Colorado 80303, USA 
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APPENDIX D 



TRANSFER FUNCTION USAGE 



The spectriim analyser when in the dBV display mode 
reads: 

N(dBV)=201og • 

0 dB referenced to 1 nT and not 1 volt is desired. Thus: 



. A volts . , X volts/nT '/ A volts . / X volts . 

^1 volt ' ^X volts/nT' ^X volts/nT' ^1 volt nT^ 

where B (the field measured by the system) is: 

R= f A volts > 

^X volts/nT^ 

so referenced to 1 nT we have: 

N(dB )=201og + 20log X . 

X is the system transfer function given in units of 
volts/nT. Since spectrum analysers use various averaging 
windows and have various numbers of bins the following 
correction must be applied: 

B^=$(BW) , 



and: 



number of bins correction 



Now we have: 



N(dB) = 201Og[$^'^^ (BW) ^'^^]+201Ogx 
and finally after rearranging: 
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101og$ = N (dB) - 201ogX - lOlog(BW) , 



where 101og$ is 
density. $ has 



the corrected reading or the power spectral 

2 

the units of nT /Hz. 
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